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Abstract. We consider the XXX homogeneous Gaudin system with N sites, both in classical 
and the quantum case. In particular we show that a suitable limiting procedure for letting 
the poles of its Lax matrix collide can be used to define new families of Liouville integrals 
(in the classical case) and new "Gaudin" algebras (in the quantum case). We will especially 
treat the case of total collisions, that gives rise to (a generalization of) the so called Bending 
flows of Kapovich and Millson. Some aspects of multi-Poisson geometry will be addressed 
(in the classical case) . We will make use of properties of "Manin matrices" to provide explicit 
generators of the Gaudin Algebras in the quantum case. 
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1 Introduction 

The Gaudin model [21] was introduced by M. Gaudin as a spin model related to the Lie algeb- 
ra 5/2; and later generalised to the case of an arbitrary semisimple Lie algebra g. 
The Hamiltonian is 



dim g 
a=l i^j 



where {xa}, a = l,...,dim3, is an orthonormal basis of g with respect to the Killing form 
(and x"" its dual). These objects are regarded as elements of the polynomial algebra S{g*)^^ in 
the classical case, and as elements of the universal enveloping algebra U(q)^^ in the quantum 
case, as 



.(0 



1 



(g) Xa 



i-th factor 

Gaudin himself found that the quadratic Hamiltonians 

dim i 
k^i a=l 



^{»)™a(fe) 



Zi 



Zk 



*This paper is a contribution to the Proceedings of the XVIIth International Colloquium on Integrable Sys- 
tems and Quantum Symmetries (June 19-22, 2008, Prague, Czech Republic). The full collection is available at 
|http://www.emis.de/journals/SIGMA/ISQS2008.html| 
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provide a set of quantities that commute with Hq, for any choice of pairwise distinct points 
zi, . . . jZtv in the complex plane. The corresponding classical quantities are "constants of the 
motion" for the Hamiltonian flow generated by Hq- For instance, in the s/2 case, we have 

N 

Hg= Y1 ^^'^^^^^ + e^'V^^'^ + f^'^e^^\ 

^ /,W/,(fc) + e«/(fc) + /WeW 
Hi= 2^ , i = l,...,N. 

Later it was shown (see, e.g., j23j ) that - in the classical case - for a general (semisimple) 
Lie algebra g the spectral invariants of the Lax matrix 

^ Z- Zi 

t,a 

encode a (basically complete) set of invariant quantities. We recall that under the name spectral 
invariants of Lq we mean the quantities 

Res(z-Zi)^"^a(-^^G), i = l,...,N, 

Z=Zi 

where is a complete set of Ad-invariant quantities for g, and the exponents run in an 
appropriate set. For instance, for g = gl{r), which is the case we will basically consider, we can 
choose ^"{Lg{z)) to be Tr (Lg{z)°^) (or, equivalently, the coefficients pa of the characteristic 
polynomial of Lg{z)), with q = 1, 2, . . . , r; in this case, ka will run from to a — 1. 

The meaning of "basically" complete refers to the fact that, e.g., the residues of the trace 
of Lq{z) at the different points Zi, i = 1, . . . , N are not independent, since Lq{z) is regular at 
infinity. Actually, the spectral invariants of Lg{z) generate the Poisson-commutative subalgebra 
of maximal possible transcendence degree in the algebra of diagonal invariants S{q* © • • • 00*)^. 
To obtain a maximal Poisson-commutative subalgebra in ^(g* © • • • © g*), one has to add 
^(dimg[„ + rankg[„) independent Hamiltonians (the Mischenko-Fomenko generators for the 
diagonal g). This is due to the fact that, being the models we are considering homogeneous, 
they are invariant with respect to the diagonal adjoint action of the group G. We will later refer 
to this fact as the global G invariance of the models. The enlarged set herewith defined gives to 
the (classical) Gaudin system the structure of a Liouville integrable system. For instance, in the 
case g = 5^2) one has to complement the residues of the trace of at its singular points zi with 

N 

a further quantity, e.g. a component of the total "spin", say Sz = Yl should however be 

j=i 

noticed that the number of integrals to be added to the spectral invariants equals, independently 
of the number N of "Gaudin generalized magnets", the rank of g. So, with a slight abuse of 
language, we will sometimes speak of a complete set of integrals referring only to the spectral 
invariants. 

The quantum case was the subject of intense investigations. It is outside of the size of 
this paper to historically frame such a research line. For our purposes, it is of paramount 
relevance the fact that, in [18], Feigin, Frenkel and Reshetikhin proved the existence of a large 
commutative subalgebra A{zi, . . . , zn) C U{q)®^ containing the quadratic elements Hi (see, 
also, [lOl [m mi [20I EOl EH). For g = s/2, the algebra A{zi, . . . , zat) is generated by the Hi, the 
additional global element Sz and the central elements of C/(g)®^. 

In the other cases, the algebra A{zi, . . . , z^) has also some new generators known as higher 
Gaudin Hamiltonians. Their explicit construction for g = gl{r) was obtained by D. Talalaev 
and further discussed in papers by A. Chervov and D. Talalaev [7l[8]. 
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In the present paper we will discuss the problem of considering - both from the classical and 
from the quantum point of view - what happens when the arbitrary points zi, . . . , ztv appearing 
in the Lax matrix, and in the quadratic Hamiltonians Hi glue together; in particular, we will 
pay special attention to the "extreme" case, when in some sense all points collide; indeed, this 
extreme case gives rise to (the generalisation of) a remarkable integrable system, called the 
"bending flows" systems introduced and studied - for the case q = so(3) - by Kapovich and 
Millson [23] in their work concerning moduli spaces of polygons in M^. 



2 The classical case 

In the classical case, to shorten notations, it is customary to write the Lax matrix Lg{z) (jl.ip 

TV 

as Lg{z) = i~F- satisfies the Poisson algebra 



=1 

n 



{Lg{z)® 1,1 ®Lg{u)} 



,Lg{z)(S)1 + 1®Lg{u) 



where 11 is the permutation matrix Ii{X ®Y) = Y®X. 

As it is well known, this linear r-matrix structure (together with suitable reductions) is 
associated with a huge number of classical integrable systems, such as Neumann type systems 
and the n-dimensional Manakov tops, finite gap reductions of the KdV equations, as well finite 
gap reductions of its generalisations to arbitrary simple Lie algebras known as Gel'fand-Dickey 
hierarchies, Hitchin's system on singular rational curves and so on and so forth (see, e.g., [291 [T] 
and the references quoted therein). 

Also, it is well known that the r-matrix Poisson brackets presented above can be seen as 
a kind of "shorthand notation" for the following situation: 

• The phase spac^ of the A^-site classical Gaudin model is g*®^. 

• The "physical" Hamiltonian is a mean field spin-spin interaction, 

1 ^ 

• The Poisson brackets defined by the r-matrix formula are just product of Lie-Poisson 
brackets on g*®^. 

N 

• The definition of the (classical) Lax matrix as Lg{z) = ^ — ^ defines an embedding of 

our phase space q®^ into a Loop space Lq ~ q{{z))- 

We now turn to the discussion of limits of the above situation when some of the points zi, . . . , z^v 
collide (or glue together) in a suitable sense. We shall discuss both algebraic aspects (that is, 
the corresponding Lax matrices) and address also some Poisson aspects. 



2.1 Lax matrices and their limits 

To consider limits of the Gaudin algebras when some of the points zi, . . . , ztv glue together, we 
can proceed as follows. We keep some (say the first k) points zi,...,Zfc "fixed", and let the 
remaining N — k points glue to a new point w, via 

Zk+i = w + sui, i = l,...,N-k, Zi^Zj, Ui^Uj, s ^ 0. (2.1) 

We can represent this procedure as in Fig. [H 

'^In this paper we will always deal with metric (or reductive) Lie algebras (and, in particular, with gl{r)), so 
that we will tacitly henceforth identify g with g*. 
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To see what happens in this hmit we first observe, that obviously 

N 

k E 

Lg{z)^L2{z)=Y,^ + '-^^^^, s^O- (2.2) 

^-^ Z — Zi z — w 

1=1 

Considering only this limit would however be too naive, and indeed, we see that (even in the 
case of g = gl{2)) the number of Hamiltonians obtained from the spectral invariants of (|2.2|) is 
not sufficient to yield complete integrability. 

One can notice that, in the gluing procedure, some residues of the Lax matrix should also be 
taken into account. That this is the case can be seen by a proper rescaling argument as follows. 

Let us introduce a new variable z = {z — w)/s, so that z = w + sz, and rewrite the Lax 
matrix in terms of the new variable as 

k „ N 



Lg{z) = \] ^ + V 



W + SZ — Zi W + sz — W — SUi 

1=1 * i=k+l * 

Multiplying this by s we can see that another Lax matrix appears for s — > 0, that is, 
/ k N \ N 

y^^^yj + sz-z, ^^_^^s{z-u,)) ,^^,Z-Ui 

The justification of this rescaling and limit is simple. Indeed, for s 7^ 0, the ring of Hamiltonians 
obtained from sLg{z) coincides with that obtained from Lq{z). Thus the Hamiltonians obtained 
from the limiting Lax matrix Li{z) are limits of the corresponding Hamiltonians of the Gaudin 
model. 

Finally we conclude that starting from the Lax matrix with generic (distinct) points zi, . . ., zat, 
we can associated, to the gluing of (|2.ip the following pair of "Lax matrices": 

JV 

Li(z)=V^ and L,{z)=y^^ + '^^^. (2.3) 

Z-Ui ^Z-Zi z-w 

Notice that, in the general case, we can choose the gluing procedure to be explicitly given by, 
e.g., 

Zk+i = w + s{zk+i-w), sG(0, 1) 
and, using invariance w.r.t. transformation of the spectral parameter z z — w, trade the 

TV 

matrix Li of (US} for Li{z) = Yl 
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In particular, in the example of the Fig. [H we would associate, to the "generic" Lax matrix 

5 

Lg(z) = Y] the two matrices 

1=1 

Li{z) = \ \ , and L2{z) = \ \ . 

Z — Z^ Z — Zj^ Z — Z^ Z — Zl Z — Z2 z — w 

Proposition 1. For every choice ofw G C the family of spectral invariants 7i^^\ Ti.^"^^ associated 
with the Lax matrices Li and L2 satisfy the following properties: 

1. The elements ofH^^\ H^^) commute w.r.t. the standard (diagonal) Poisson brackets on q . 

2. The dimension of the Poisson commutative suhalgehra Tti^2,w generated by the spectral 
invariants ■H(^) and ■^(2) ^hich are those, respectively, associated with the Lax matri- 
ces Li{z) and L2{z), coincides with that of the spectral invariants associated with the 
generic Lax matrix. 

3. The physical Hamiltonian Hq lies in 7ii^2,w 

Proof. The commutativity property (the first assertion) holds thanks to the following facts. 
First, it is clear that if hi, /12 are taken either both from the subfamily T^*^^-* or the subfamily 7^^^^ 
the assertion is trivially true. Let us suppose thus that hi S TY*-^-* and /12 S TL^"^^ . The fact that 

{/ii,/i2} = follows from the fact that /12 depends on Xi,X2, ■ ■ ■ ,Xk only through the sum 

k 

Xfc = ^ Xi, and hi is invariant under the diagonal action of g. 
j=i 

The second assertion holds thanks to the results of Jurco briefly mentioned before about the 
number of functionally independent spectral invariants of a Lax matrix of Gaudin type, as well 
as to the functional dependence of L2{w) on the "variables" Xi, . . . , X^- 

The third one can be checked by induction. ■ 

Remarks. 1) The meaning of this proposition is that, once the global G-invariance has been 
taken into account, the procedure of glueing discussed above provides new families of Liouville 
integrals for the "physical" Gaudin Hamiltonian II g- 

2) It is fair to say that a procedure somewhat similar to this was presented by the late 
V. Kuznetsov in [25]. However, possibly owing to the fact that the author restricted himself to 
the case of rank 1 (namely, mostly to the case g = so{2, 1)), he did not consider the need to add 
the further singular point w. This is crucial to obtain, in the general case, the correct number 
of Liouville integrals for the "new" system. 

3) Our glueing procedure differs substantially from the one considered in [27] (see also [2]), 
where the limiting procedure produces higher order poles in the Lax matrix, and is associated 
with non-semisimple Lie-Poisson algebras, that are Inonii-Wigner contractions of the original 
one (see, for our case. Section [2^ . 

2.2 On the Poisson geometry of the hmiting procedure 

The interplay between classical Lax matrices. Loop algebras and Poisson manifolds is nowadays 
well known, and was settled mainly in the works of the Leningrad's school (see, e.g. the review 
by Reyman and Semenov-Tyan-Shanski in |29j). It is encoded in the notion of i?-operator as 
follows. 

On the space q*{{z)) of Laurent polynomials with values in (the dual of) a Lie algebra g, 
there is a family of mutually compatible Poisson brackets, {■,-}k associated with a family of 
classical /^-operators 



R,{X{z)) = {z^X{z))_^-{z'X{z))_ 
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via the formula 

{F,Gh{X) = {X, [i?,(V(F)),V(G)]) - {X, mV{G)),V{F)]). 

Spectral invariants of a Lax matrix (say, polynomially dependent on the spectral parameter z) 
form an Abelian Poisson subalgebra w.r.t. any of the brackets {•, ■}k- 

Let us consider the space of A^-th order polynomials q*j^ = 5*[[z]] mod z^~^^. It can be shown 
that (see, e.g.. 



The brackets associated with Rq, . . . ,i?Ar on q*{{z)) restrict to the affine subspace 
Q%A--=lxeg*{{z))\X{z)=z''- 



N >| 

-^A + ^Z^XA, 
1=0 J 



where ^ is a fixed element of q, and thus on q*^ = q*^q, that is when ^ = which is the 
case we plan to consider in this paper. 

• These brackets are mutually compatible and give rise to multi-Hamiltonian structures on 
the finite dimensional manifolds flj^r^- A straightforward observation [T^ is the following: 

N 

we can associate with any polynomial of degree Qiz) = ^ Poisson bracket 

i=0 

(of the RSTS family) via: 

TV 
i=0 

In this setting, we can recover the standard (product) Lie-Poisson structure on g® ^ as follows. 
Considering the "Lax m 
framework, by the sum 



Considering the "Lax map" Lq = Y] —z- the diagonal structure can be obtained, in the RSTS 

1=1 



N-1 

{•,-}5 = {-,-}iV+E(-l)^^^{-'-}- 

j=0 

the (Tj being the elementary symmetric polynomials in the quantities zi, . . . , zn. That is, the 
standard r-matrix Poisson bracket can be regarded as the one associated with the polynomial 

N N-1 

S{z) = lliz - z,) = + ^ (-l)V,z^-i. 

1=1 i=0 

As a side remark, one can also notice that this gives the possibility of constructing a bi- 
Hamiltonian structure for the classical Gaudin case. Indeed the structure defined by the poly- 
nomial (^^Y^^ provides, together with the standard one, a bi-Hamiltonian structure for the 

Gaudin model, such that the spectral invariants of Lg{z) fill in recursion relations of (genera- 
lized) Lenard-Magri type. More in detail, suitable combinations thereof give rise to anchored 
Lenard-Magri sequences, in the terminology of Gel'fand-Zakharevich |22j . 

What is more important for the present paper is that this picture suggests and gives the 
opportunity of studying, from this point of view, the limits of (suitable combinations of) the 
RSTS Poisson structures when some (and, iteratively, possibly all) poles of the Lax matrix 
zi, . . . , 2:jv glue together. 
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This means the following: if we pull back, via the map 

N 



n( 

i=l 




the RSTS Poisson brackets {•, - jfc, k = 0,. . . ,N, we obtain families of linear Poisson brackets 
on Q*j^ = Q*[[z]] mod z^~^^ that depend rationally on (zi, . . . , z^). 

In this framework, a natural question arises, that is, whether one can associate suitable 
Poisson structures to the procedure of gluing (some of the) points Zi discussed - in the Lax 
setting - in Section 12.11 In the example of the limit pattern depicted in Fig. [H one of the 
possible Poisson structure^ that can be obtained is represented by the operator (for the sake of 
simplicity, we set w = Z3, g = gl{r)) 



V 





223 [-^1, 







223 [^1,-] 
[^22,-] 

-212 [^3, 

-212 [^4, 
-212 [^5, 





-212 [X3 



213234 





-212 [^4 
213234 





-212 [^5, 



245 
213235 

^45 



[^3 
[^3, 



245 
[-P44, 



[^3 



213235 
245 

[A5, 



[^3 



45, 



[Pi 



55, 



(2.4) 



where 



-P22 
P44 



Zj, and 



Z23 {X2 - Xi) + Z12 (X3 + X4 + X5) 



213234 

245 
213 



235 — 245 
245 



213234 
245 



X, 



P45 = 4^ (215X4 + 214X5), 



'45 



55 



£13235 
245 



Y _L 35 V _L ^34 - 245 

A3 H A4 H 

245 245 



X, 



Using the fact that both tangent vectors and one-forms on the phase space of this five-site model 
can be identified with five-tuples of matrices the meaning of the representation (j2.4p (see also 
the subsequent equations (j2.5p . (|2.6p can be illustrated as follows. 

The Hamiltonian vector field associated via V to the one-form (ai, . . . , 05) is given by Xi = 

5 

^ [Pij^aj]. E.g., Xi = Z23[Xi,Q;2], and so on and so forth. 

3=1 

It is not difficult to check that the spectral invariants of the matrices 



X3 



23 



+ 



X4 



24 



+ 



X. 



and 



24 



Xi 

; - zi 



+ 



X2 



+ 



Xi + X2 + X3 



22 2 - Z3 

(besides, as proven in 



do commute with respect to the Poisson brackets V defined by ([2 
Proposition [H commuting w.r.t. the standard Poisson brackets). 

Conjecture. For every limit pattern vr described in Section [2?T] there are suitable linear combina- 
tions {•, 'Itt of the Poisson brackets above that remain regular. Endowing the phase space g*®^ 
with these limiting brackets, one obtains bi-Hamitonian (or, possibly, multi-Hamiltonian) mani- 
folds. The integrable systems defined, (according to the Gel'fand-Zakharevich scheme) by these 



^There is two-parameter family of such structures. 
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multi-Hamiltonian structures admit Lax representations, whose Lax matrices are exactly those 
constructed in Section [2.11 In other words, the spectral invariants of the Lax matrices associated 
with the limiting pattern vr patterns provide, w.r.t. the Poisson pencil formed out of the diagonal 
bracket and the bracket {•,-}n Lenard-Magri sequences (possibly in a generalised sense). Also, 
the brackets {•, -Itt are semisimple Lie-Poisson brackets. 

This conjecture has not been fully proven yet, but we have checked it in a significant number 
of examples. We notice that, in the "extreme gluing case", that is, when all Zi glue according 
to the pattern Z2 — > zi, z^ Z\ (= Z2), 2^4 z\ (= z^ = z^), ... we can obtain the following 
"limit" bracket which is independent of any parameter: 

{F, G}ii,nit = UjkTr {VFi [Xk , VG, ] ) with 

rijk = {k — l)6ij6jk — 0(i_j^'^5ij + 9(^j-i)5ik + 6[i-j)5jk-i 

where VFj, (resp. VGj) represents the differential of F (resp. G) w.r.t. the i-th entry Xj. 

As it has been shown in [13] these Poisson brackets are still compatible with the standard 
r-matrix ones. 

To elucidate the formula above, we notice that, in the 4-site case we have the following 
representation - by Poisson operators - of the brackets as follows: 



/ [^1 



-^standard r 



limit 



V 



V 



















[X2-Xi,-] 




[^2,-] 






[2X3-X2-Xi,-] 


[x,A 






[^3,-] 


[3x4 - E ^> 

i=l 



J 



(2.5) 



(2.6) 



According to the Gel'fand-Zakharevich (or Lenard-Magri) scheme, Pstandard r and Pumit define 
the gl{r) (as well as g) - valued generalisation of the so called Bending flows, introduced in the 
case corresponding to g = 0[(2) by Kapovich and Millson [25, and further discussed in |16j . 



2.3 The Bending system 

Bending flows are defined on the moduli space Aiy. of (A^)-gons with fixed sides lengths r = 
(ri, . . . ,r/vr), > 0. In [2?] it was shown, among other properties, that: 

• Mr is a smooth (2A^ — 6)-dimensional open manifold, whose compactification is achieved 
when some of the side lengths vanish. 

• Mr is endowed with a natural symplectic structure, since it is a symplectic quotient of 
products of spheres: 

{(ei, 62, ... , cn) e Sl^ xS^x-.-x S^.J//SO{3), 

where the symbol // denotes the symplectic quotient on the null manifold of the moment 
map defined by the diagonal S0{3) action on the spheres Sr^. 

On this moduli space of polygons a completely integrable Hamiltonian system can be defined, 
indeed called the Bending system or system of bending flows. 
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Figure 2. Polygon in M"^. Figure 3. The Bending flows for polygons ((so(3))-case). 

The action variables for such system are the lengths of diagonals stemming from one vertex, 
and the angle variables are (indeed) the dihedral angles (see Figs. [2] and E]). Geometrically, the 
flows bend one part of the polygon around the corresponding diagonal, keeping the rest fixed, 
whence the denomination. 

In the papers [HI [13] the following picture was established: g-Bending flows can be defined 
on the same phase space of the g-Gaudin system, the link between the two (in the Kapovich 
Millson A^-gon case) being the fact that is a symplectic leaf of so(3)^, and the condition for 
the polygon to close is just the choice of the special value = of the image of the momentum 
map. Also, the following results were established: 

• g-Bending flows admit a set of — 1 Lax matrices of the form 

N 

Lkiz) = zXk + J2 k = h---,N-l. (2.7) 

i=fc+l 

• Hamiltonians come in "clusters" , each cluster being associated with the corresponding Lax 
matrix. For g = gl{r) 

Hk,m = l^l^T^Lf{z). (2.8) 

• Separation of variables can be performed in this scheme. Separation canonical conju- 
gated variables come in "clusters", associated, via the bi-Hamiltonian version [15] of the 
Sklyanin magic recipe [33], to each of the Lax matrices. Actually, it has also been shown 
that integration of the Hamilton-Jacobi equation involves Abelian differentials of order 
independent of the number of sites. Moreover, for g = sl{r) the genus of the underlying 
Riemann surface can be computed to he g = , 

For further use, we remark that the ring of spectral invariants of a Lax matrix is, as it is well 
known, invariant w.r.t. rational changes in the spectral parameter and/or multiplication by 
a rational function of the spectral parameter. Thus, instead of the Hamiltonians (j2.8p asso- 
ciated with the "parameter free" Lax matrices (j2.7p . we can consider, as an equivalent set of 
Hamiltonians for the generalised gl{r) bending flows, the spectral invariants associated with the 
matrices 

X2 Xi X^ Xi + X2 

J^l — 1 , — 1 , ■ ■ ■ , 

Z — Z2 Z — Z\ Z — Z2 Z — Zi 
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Figure 4. Triangulations and trees. 



Af-l 

L^-i = ^ + ^^, (2.9) 

Z — Z2 Z — Z\ 

rationally dependent on z. We finally want to add the following observation. As it is pictorially 
suggested in Fig. U there is a natural interpretation of our procedure in terms of triangulations 
of a polygon and graphs. 

Actually, we have already seen (see Fig. [1]) that we can interpret our procedure of glueing 
points in terms of building up tree-graphs out of the set of singular points of the generic Lax 
matrix of the Gaudin model. In particular, the glueing pattern depicted on the right of Fig. H] 
corresponds to the limit in which we produce the "Bending" integrals, and, in the topological 
point of view, exactly gives (provided we identify the point z\ with the root), the dual graph to 
the triangulation of the right hand side of the picture. 

It is outside of the size of the paper to fully discuss these issues, which however the authors 
think being of a certain interest. It should also be noticecH that tree diagrams corresponding to 
the glueing patterns described in this paper appear in the theory of moduli of pointed rational 
curves. 



3 The quantum case 

3.1 FFR Gaudin algebras and their limits 

The existence of a large quantum commutative subalgebra A{zi, . . . , zat) C U (g)®^ containing 
the Gaudin Hamiltonians Hi was shown in [18] . As we recalled in the introduction, whenever g 
has rank greater than one, A{zi, . . . , zn) has new generators, besides the quadratic ones, known 
as higher Gaudin Hamiltonians. 

Roughly speaking, the definition of A{zi, . . . , zn) can be obtained identifying a commu- 
tative subalgebra of the enveloping algebra U{q t~^C[t~^]) as follows. To any collection 
zi,. . . ,Zn of pairwise distinct complex numbers, one can naturally assign the evaluation map 
U{q (S) t~^C[t~^]) — > [/(g)®^. The image of the commutative subalgebra in U{t-^Q[t-^]) under 
the composition of the above homomorphisms, denoted by A{zi, . . . , Zn), is called (quantum) 
Gaudin algebra. 

It is useful to recall the following well known facts from the theory of Lie algebras. Let g be 
a Lie algebra. The universal enveloping algebra U{q) bears a natural filtration by the degree 
with respect to the generators. The associated graded algebra is the symmetric algebra S{g) by 
the Poincare-Birkhoff-Witt theorem. To every element ^ in [/ (g) there corresponds uniquely its 



^We thank G. Felder for an observation concerning this point. 
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image ^ in S{q), since S{g) corresponds to the graded algebra associated with U (g). This element 
is customarily called, in the theory of Gaudin and Heisenberg spin systems/chains, the "classical 
limit" of ^. Also, whenever g is reductive, S{q) ~ S(q*), and U(g e • • • g) = (t/fg)*^^). 

V ' 

A' times 

The problem of finding explicit representatives for generators of A{zi, . . . , zm) in the case 
g = gl{r) was concretely solved by Talalaev some ten years later after the paper by Feigin, 
Frenkel and Reshetikhin [TS], by means of the following construction. 

Theorem 1 (Talalaev, 2004). Let L{z) be the Lax matrix of the gl{r)-Gaudin model, that is, 
let L{z) satisfy the r-matrix commutation relations 



[L{z) ® 1,1® L{u)] 



^ -,L{z) 01 + l®L{u) 



u 



(3.1) 



Consider the differential operator in the variable z "DET" (9^ — L(2)) = ^ QHi{z)d\; Then: 

i=0,...,r 

Vi,j € 0, . . . ,r, and u,vGC, [QHi{z)\;,=u,QHj{z)\:,=^] = 0. 

The QHi{z)'s generators of (a full set of) quantum mutually commuting quantities. We shall 
discuss the meaning of the symbol "DET" appearing in the formulation of Talalaev's theorem 
in Section 13.21 where we shall also comment on its proof. 

For the moment, let us remark that the Gaudin algebra(s) defined by Feigin, Frenkel and 
Reshetikhin - and concretely identified by Talalaev's formula for A'^-spin gl,, Gaudin systems - 
explicitly depends, in general (that is when A > 3) on the points zi, . . . , zn- In another parlance, 
its has moduli. However it holds (see for the proof) the following 

Proposition 2. A{zi, . . . , zjy) is invariant under permutations and under simultaneous resca- 
lings Zi azi + [3, a,(3 (z C; thus the "two site" algebra A{zi,Z2) is independent of zi, Z2. 

Let us now consider limits of the quantum Gaudin algebras when some of the points zi, . . . ,zm 
glue together according to some pattern. Here, obviously enough, we will follow the same glueing 
procedure as in Section [2l 

We recall that the basic procedure is that we keep some (say the first k) points zi, . . . ,Zk 
"fixed", and let the remaining N — k points glue to a new point w, via 

Zk+i = w + sui, i = 1, . . . ,N - k, Zi^ Zj, Ui^ Uj, s ^ 0. 
To describe this limit in the quantum case, it is useful to introduce the maps 

Dk,N ■■= id^' ® diag^v-fc : C/(g)«(^+^) ^ C/(g)«^, 
defined by 

Dk,N{Xi ■ ■ ■ (g) Afc (g) Afc+i) = Ai • • • Afc (g) Afc+i (g • • • Afc+i 

^ V ' 

N—k times 

and 

Ik,N ■■= 1®^ ® id^(^-^) : i7(g)®(^-'=) ^ Uig)®^, 
defined by 

Ik,N{Xi g) •••AAT.fc) = 1 (g) ■ ■ ■ (g 1 g)Ai (g) ■ ■ ■ (g) XN-k- 

k times 

Furthermore, let us define the algebra 

A(^z^^,„^Zk,w),{ui,...,UM-k) '■= Dk,NiA{zi, ...,Zk, W)) ■ Ik,NiA{ui, UN-k)) 

In [1] it is proven the following 
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Theorem 2. The algebra ^(2j,...,zj.,u,),(ui,...,njv-fc) commutative; 

limA{zi, ...,Zk,W + SUi, ...,Z + SUN^k) = -^{zi,...,Zk,z),{ui,...,u,^^i:)- 

Remarks. 1) In some sense we arrive at a kind of factorisation of the Umit algebra by "adding" 
one point. Indeed, the commutative algebra ■A.(^zi,...,Zk,w),{ui,...,UM-k) evolves the FFR algebra as- 
sociated with the points {zi, . . . , Zk,w), and a FFR algebra related with the points ui, . . . ,uj^-k- 
Iterating this limiting procedure described above we can obtain some new commutative subal- 
gebras in U{q)^^ , just like in the classical case we found Poisson commutative subalgebras 
of 5(0)®"^. 

2) In the case g = gl{r) the classical limits (in the sense specified above) of, respectively, 
7fc^Ar(^(iii, . . . , UN-k)) and Di^^n{A{zi, . . . , Zk,w)) coincide with the ring of the spectral invari- 
ants of the matrices 

TV 

Li = and L2{w) = h — 

^ — ' r. — II.- ^ — ' 



Z — Ui ^ — ' Z — Zi Z — W 

i=k+l 1=1 



already discussed in Section [2l 

In particular, as in the bending flow case, we can iterate the procedure to pass from A{zi, Z2, 
...,zn) to A(zi,z2,...,ZN-2,w) (zjv_i,2jv), and finally obtain the subalgebra 

Aim ^ C Uisr^. (3.2) 

In the process, we have to use translation invariance of A{zi, . . . , zm), as well as the property 
that the two site Gaudin algebra A{u,v) is independent of {u,v). 
The limit algebra (j3.2p is generated by 

Di,NiAizi,Z2)), l®L>l,7V_l(^(zi,Z2)), 1^(^-2) 0^(zi,Z2). 

More explicitly: 

Proposition 3. The subalgebra Aymi is generated by elements Hi'^ € U {gl{r))®^ such that their 

— (a) 

classical "limits " , are given by 



i?iJ(Xi,...,XO :=-Res-^Tt[Xk + z ^ X, 





k+l / 

In other words, the classical limits of the (still unspecified) "quantum Bending Hamiltonians" 

N 

do coincide with the spectral invariants of the Lax matrices Li.{z) = zXk+ ^ Xj, i = 2, . . . , N , 

associated with the bending flows. In turn, these can be traded for the spectral invariants of 
their "rational" analogs 

Li{z) = \ , L2{z) = \ , 

Z — Z2 Z — Z\ Z — Z2 Z — Z\ 

N 

Lm-1 = + (3.3) 

Z — Z2 Z — Z\ 

In the last part of the paper we shall show how, using Talalaev's results briefly reminded 
before as well as some of the results of j3] we can find a suitable quantization of the traces 
of the powers of the Lax matrices (j3.3p (as well as of, more generally, of the quantum Lax 
matrices L\{z\ L2{z) associated with an elementary "gluing" procedure. 

To this end, we need to introduce the notion of Manin matrix. 
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3.2 Manin matrices and the quantization of traces of Lax matrices 

We consider the (differential operator valued) matrix considered in Talalaev's Theorem [H that 
is, M[z) = dz — L{z). The standard linear r-matrix commutation relations (j3.ip imply that the 
matrix elements of M. satisfy special commutation relations. These are the defining commuta- 
tion relations of a class of matrices (with non-commutative entries) called in [3] Manin matrices. 
The terminology originates from a seminal paper of Yu.I. Manin on quantum groups [26]. 

Manin matrices are, in a suitable sense, matrices associated with linear maps between com- 
mutative rings. Their (more operative) definition can be given as follows: 

Definition 1. Let Mij be a matrix with elements in a non commutative (unital) ring TZ] we 
call it a (column) Manin matrix if: 

• elements in the same column commute among themselves; 

• commutators of the cross terms in any 2x2 submatrix are equal: 

[Mij ,Mki] = [Mkj , Mil] , e.g. [Mil , M22] = [M21 , M12] and so on and so forth. 

Manin matrices admit a natural definition of determinant. Indeed, if M is any matrix, one 
can in principle define "a" determinant of M by column expansion: 

rv 

det M = det^°i M=Y,i-^r H ^-(0,- 

trGSn j=l,...,n 

where Sn is the group of permutations of n letters, and the symbol r\ means that in the product 
Y\ M^(^i^^i one writes at first the elements from the first column, then from the second column 

i=l,...,n 

and so on and so forth. 

In the case of Manin matrices the characteristic property of total antisymmetry of the deter- 
minant is preserved: 

Proposition 4. The column determinant of a Manin matrix does not depend on the order of 
the columns in the column expansion, i.e., 

ypeSn det-iM= ^^(-1)'^ H M,(p(,)),p(,). 

crGSn i=l,...,n 

In particular, a good notion of "DET" in the formulation of Talalaev's quantum integrals of 
the motion is given by the column determinant of M{z). 

Among the properties of "ordinary" linear algebra that Manin matrices satisfy, we mention 
the following: 

• The inverse of a Manin matrix M, whenever it exists, is still a Manin matrix. 

• The (left) Cramer formula holds: 

M^m = det™'(M)l, 

where M^'^^ is computed - via column determinants - in the same way as in the ordinary 
case. 

• Schur's formula for the determinant of block matrices holds: 

det™' ( ^ ^ ) = det^°\A) det^°\D - CA-^B) = det^°\D) det™' {A - BD-^C). 
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• The Cayley-Hamilton theorem, det™ (t — M)\f-M = 0, holds true. 

Further properties of Manin matrices, as well as full proofs of these facts can be found in [5]. 
Particularly important, for our purposes, is the the fact that a suitable redefinition of the traces 
of the powers of the quantum Lax matrix allows us to trade coefficients of the characteristic 
polynomials for these "normalised traces" . The key property concerns the quantum counterparts 
of the classical Newton identities between these two classes of Ad- invariant quantities. The 
coefficients ai of the characteristic polynomial of a n x n matrix M are, as it is well known, the 
elementary symmetric polynomials in the eigenvalues Aj of M; however in the ring of symmetric 
polynomials in the n variables Aj, one can consider other sets of generators and in particular the 
so called power sums 

n 

ri = Y,Xi = 'n^{M'), i = l,...,n. 

k=l 

In the case of matrices with commuting entries, the family {cj}, i = l,...,n and {rj}, i = 
1, . . . , n are related by the identities (the Newton identities): 

{-l)''+^kak= Yl i-^y^^n-i■ (3-4) 

i=0,...,fc-l 

The same identities hold true, in the case of Manin matrices, for the corresponding quantities. 

Theorem 3. Newton identities of the form (13. 4p between Tj = TrM'^' and the coefficients ai of 
the expansion o/ det'^°'(t + M) in powers oft hold. 

Proof. First one can observe that the following property, whose proof is straightforward, 

Tr {t + M)'^^ = dt det^°\t + M) 
hold. Then we have that 

i ^ Tr((-AfA)'=)=Tr^ = Tr((det-'(t + M))"^(t + A/r^j) 

fc=0,...,oo 

= (det^°'(t + M))"^Tr(t + M)^^j = ( det'=°'(t + M)) "^9* det^°'(t + M). 
The result is obtained substituting —M for M in these formulas. ■ 



The quantization of the traces of powers of the Lax matrices can be defined according to the 
Newton identities and the following considerations. First, a straightforward computation that 
makes use of the the r-matrix commutation relations (13. ip . proves the following 

Proposition 5. Lei0 

'-^T^ ,3.5) 

be the Lax matrix for the quantum Gaudin model. For any number of sites N , and different zi 
the matrix dz — L{z) is a Manin matrix. 



''in these and subsequent formula we use the symbol "to stress the fact that in the quantum case Xi are to be 
regarded as elements in the universal enveloping algebra of gl{r). 
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It might be mentioned that this is not a specific property for the Gaudin system herewith 
considered. Indeed, Manin matrices enter other topics in (quantum) integrability; for instance, 
e" "Tgi(r)-Yangian(^) IS Manin, where 

-^g;(r)-Yangian(^) the Lax (or "transfer") matrix for the 
Yangian algebra Y{gl{n)), satisfying quadratic r-matrix commutation relations (see [3] for more 
details). 

The "quantization" of the traces of a Lax matrix of the form (13. 5p can be achieved as fol- 
lows. From Talalaev's theorem we know that the coefficients QH[z) of characteristic polynomial 
det'^°^(5^ — L{z)) do commute (at the quantum level), and from the theorem about Newton iden- 
tities we know that we can trade these coefficients for traces of the powers of the corresponding 
Manin matrix. 

Basically, we simply have to remark that we should not consider the quantities Ti{L{z)^)^ 
which, as it has been shown in [6] do not, generally speaking, commute, but rather the traces 
of powers of the Manin matrix 

k 

TV((5, - L{z)f) = Y,{Q T^)]{z)dl-^, k = l,...,r. (3.6) 
i=o 

As it is easily seen, there is a recursion relation of the form Q Tr^^„]^(z) ~ Q Tr^^^(2;), and hence, 
to obtain the expected number of independent quantities, we can consider simply the coefficients 
QTr^, that is the coefficients of zeroth order of each differential "polynomial" (|3.6p . 

In turn, one sees that these quantities are given by the traces of matrices Lfc(z)["], that can 
be called "quantum powers" of L{z)] these powers are defined by the Faa di Bruno formula (see, 

e.g., m) 

ifiz) = Id, Lf{z) = Lt'\z)L,{z) -§;{Lt'\z)). 

These arguments work for any Lax matrix of "Gaudin type" , that is, of the form (|3.5p : in par- 
ticular, they hold for the (rational) Lax matrices Li{z) (j2.9p associated with the Bending flows, 
since the dependence on the spectral parameter is chosen in such a way to fulfil the canonical 
r-matrix commutation relations (j3.ip . Thus computing the traces of the quantum powers of 
each of the (quantum counterparts of the) matrices Lj(z) of (j2.9p one gets a quantization of the 
Hamiltonians of the Bending flows. 

4 Concluding remarks and open problems 

In this paper we have discussed, both in the classical and in the quantum settings, some features 
in the theory of (homogeneous) Gaudin models, and concentrated on the Gaudin algebras of 
mutually commuting quantities. We have addressed, in particular, the problem of a suitable 
definition of a complete set of Hamiltonians in the case when the poles of the Lax matrix Lg{z) 
of the model glue together. In the classical case, we have shown that a complete set of mutually 
commuting Liouville integrals can be obtained by a suitable analysis of the limiting procedure 
on Lq{z). We have then addressed some features of the (multi)-IIamiltonian properties of the 
limits of these (classical) Gaudin systems, in the framework of the Hamiltonian geometry of 
Loop algebras. Also, we touched the problem of framing the geometry of these systems within 
the combinatorial approach to the geometry of moduli spaces of pointed rational curves. In these 
respects, the results herewith collected are somewhat sketchy, and will be developed elsewhere. 

In the quantum framework, we have followed the line of [18] as well as previous results of 
ours [HE]. We reviewed how our gluing procedure for the poles of the Lax matrix of the Gaudin 
model can be applied in the quantum case to yield new quantum Gaudin algebras. Also, by 
means of results of ^ about the theory of the so-called Manin matrices, we have shown quite 
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explicitly the procedure to define, in the quantum Gaudin algebras, a set of generators that 
reduce in the classical case to the usual traces of powers of the Lax matrix. 

Besides the full solution of the problems in the Hamiltonian and combinatorial aspects of 
these systems that were just addressed in this paper, we can envisage at least two domains 
where the study of the limits discussed in the present paper might be of some interest. The 
first is the realm of Heisenberg-like models, and the theory of Yangian algebras. The second 
is the study of non- homogeneous Gaudin models, that, besides of their own interest in the 
mathematical theory of integrable (quantum) systems, proved to be important in the theory 
of strongly correlated electron systems (see, e.g., [32]). Due to the fact that non homogeneous 
systems lack global G-invariance, however, the results of this paper cannot be applied to that 
case in a straightforward manner. Work in this direction is in progress. 
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